A numerical study of the Random Dimerized XX spin-1/2 chain 



a^ 

> 
o 

in 



I 

CZ3 



I 

o 



o 
o 



X 



p. Hcnclius and S.M. Girvin 
Department of Physics, Indiana University, Bloomington, IN 47405 
(February 1, 2008) 

The effects of randomness and dimerization on the spin-1/2 XX chain are studied by a mapping to 
free fermions. Results are presented for the transverse and longitudinal components of the average 
and typical spin and string correlation functions. Contrary to previous numerical evidence, the 
decay exponents of all the above correlation functions are found to be in good agreement with the 
theoretical RG predictions for the random singlet phase. We also present cross-over functions for 
the above correlation functions in the random dimer phase. Disorder averages have been taken for 
system sizes up to N=1024, about ten times larger than in previous studies. 



I. INTRODUCTION 

It is the purpose of this work to study the effects of 
disorder and dimerization on the jgijound state of the-JCX 
spin-1/2 chain. Much analyticaltrB and nunicricalEJo 
work has been done on randonji-spiii chains, as well as 
higher dimensional spin modelslijiiZI and to begin with 
we will review some of the knowledge we have of these 
systems. |— ■ 

In 1983 Haldane pointed outEJ that integer spin chains 
have an energy gap, while half integer spin chains are 
gapless. If we introduce alternating bond strength in a 
spin-1/2 chain, anii-thereby enforce dimerization, an en- 
ergy gap is inducedlj. In the absence of disorder the spia.- 
1 and dimerized spin-1/2 chains are in the same phaseEj, 
since one can be continuously transformed into the other. 
This Haldane phase is characterized by an excitation gap, 
exponential decay of the spin correlation function and a 
non-vanishing string order parameter, to be defined be- 

Fisheiihi has recently made considerable progress to- 
wards understanding the strongly disordered random sin- 
glet (RS) phase of the spin chain. Through an asymptoti- 
cally exact RG scheme closely related to that of Dasgupta 
and Maj23 he showed that the RS phase is characterized 
by a dynamic critical exponent z — oo, meaning that 
the characteristic time scale is not a power of the chara- 
teristic length scale, but rather an exponential function, 
giving so-called activated dynamics. The RS phase has 
no energy gap and the spin correlation function between 
two given sites vary over several orders of magnitude from 
one disorder configuration to the next. The average spa- 
tial decay of the spin cortjslation functions are, however, 
described by power lawaa. Due to the broad distribu- 
tions, the typical correlations (the average of the loga- 
rithm of the correlation function) behave very differently 
from the average correlation functions. No string order 
exists. In the RS phase all the spins pair up and form 
singlets over arbitrarily large distances. Average corre- 
lation functions are dominated by these strong pair-wise 
coupling of spins, and due to the singlet nature of the 
pairing, all components of the correlation functions are 
predicted to decay with the same exponent, even if the 
underlying Hamiltonian is not rotationally invariant. If 
a weak dimerization is now enforced the system is driven 



to a random dimer (RD) phase which is gapless, like t 
RS phase, but which has a non-vanishing string orde: 
like the Haldane phase. _ 

Recently Young and RiegenlJ were able to numerically 
verify many of these striking RG predictions by using 
a mapping to free fermions. In this paper we have ex- 
tended their results to larger system sizes and we have 
studied additional correlation functions in the RS and 
RD phases. Different components of the correlation func- 
tions are found to decay with the same exponents, in 
agreement with theoretical RG predictions, but contrary 
to earlier numerical workllil. In the RD phase we present 
the full finite size scaling cross-over functions for the av- 
erage and typical values of both components of the spin 
and string correlation functions. 



II. MODEL 

We consider the following Hamiltonian H of the XX 
chain: 



N 



(1) 



where S" are spin-1/2 operators obeying periodic bound- 
ary conditions {Sn+i = Si) and Ji are positive coupling 
constants. 

In the disordered system we present results for the fol- 
lowing flat bond distributions: 



p(Jo) = {2 + d)-^e{2 + d-Jo)eiJo), 



P{Je 



(2 - d)-^e{2 



Je)0{Je). 
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where Jo and Je denote odd and even bonds. The dimer- 
ization is defined as 



[In Jo] av - [In Je]av 
var[ln Jo] + var[ln Je 



(3) 



where var denotes variance. In terms of the above distri- 
butions this reduces to 



^=iln(^ 
2 4-d 



(4) 
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In the disorder free dimerized system Jo = \ — d/2 and 
Je = l+d/2 and we use Eq. (H) to define the dimerization 
5. 

The longitudinal component of the spin correlation 
function, C^(r), is defined as 



N 
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(5) 



The transverse component, C^(r), is defined analogously. 

The longitudinal component of the string correlation 
function, 0^(r), is defined as 



N 



i=l 

+5f_^^_]_)]5f_^,.). (6) 



Using the identity 5^ — exp(z7rS'^)/2i this can be rewrit- 
ten, for spin-1/2 operators and r odd as 



N 



N 



i+r / 
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The behavior of the string correlation function in a 
dimerized system depends on whether the site index i 
in the sum above is odd or even. Therefore we have sep- 
arated the above expression into a sum over even sites, 
defining 0^{r) as 



yr+2 



N 



and a sum over odd sites, defining O^ir) as 

2r+2 _ 



Ol{r) 



N 



i odd 



{SfS^j^i ■ ■ ■ S. 



(8) 



(9) 



The difference between the the expressions is denoted 
^di''') — Ol{r) — Ol{r). The transverse component of 
the string order function is defined analoftQualy— 

This correlation function was introducec£rB'E3 to 
measure hidden long-range correlations in integer spin 
chains where the ordinary spin-spin correlation function 
vanishes exponentially. If we, for example have a per- 
fectly dimerized spin chain with J2i = oo and J2i+i — 0, 
then the spins form singlets around the strong bonds, and 
neighboring singlets are uncorrelated, with — — I 
and Of = 0. If, on the other hand there is no dimeriza- 
tion in the system — = and hence 0^ = 0. 
The string order parameter Ofj^ is defined as 



Of,, = hm Of (r) = hm Of(r). 



(10) 



We are, however, working with finite size systems with 
periodic boundary conditions and all results are calcu- 
lated at the largest distance around the ring. Results for 



C^{r) are calculated at r = N/2. C^{r) vanishes for r 
even (see Appendix) and so results are for r = N/2 — 1. 
Both components of the string order function are defined 
for odd r, and they are calculated at r = N/2 — 1, using 
our definition of the correlation functions. For simplic- 
ity we have, however, plotted all results as a function of 
n = N/2. This will introduce some corrections to scaling 
for and will be discussed below. 

All decay exponents 6 are defined according to 



f{r) cx 



(11) 



III. MAPPING TO FREE FERMIONS 

Using the Jordan- Wigner transformation we can map 
the XX model onto free fermionaHJ. Therefore we only 
need to solve the one-body problem which involves diago- 
nalizing a {N x N) matrix and hence we can study fairly 
large system sizes, up to about N = 2048 for a clean 
system, and up to = 1024 when disorder averages are 
taken over several thousand configurations. 

We apply the Jordan- Wigner transformations 



= exp(-i7r^c]cj)c 



i-l 



(12) 



St = clexp(i7r^c];cj) 
to our Hamiltonian and obtain 



H = ^ ^[cjc+i -I- Qc|+i]. 



(13) 



Next we must calculate the string correlation function in 
the fermion language: 



-J 2^ 



(14) 



3=t 



We proceed by rewriting thisE2rE3 in the form 

i+r ^ i-\-r 



(15) 



where Aj — (c] -f Cj) and Bj — (cj — Cj). This expression 
is easily evaluated using Wick's theorem. First we can 
show that 
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{AiAj) — 6ij 
{BiBj) — ~Sij 

{B,Aj) = {B,A,) = -5^, + 2(c]c, 
{A,Bi) - ~{B,Aj) = - 2(c]c,) 



(16) 
(17) 
(18) 
(19) 



and the only non-zero contractions to appear are then of 
the form {AjBi) or (BjAi). This leads to the following 
determinant: 



{A,B,) {A,B,+^) . 
{Ai+iBij (Ai+i_B,:+i) . 



^ 2 ' 

{A.,B,+r) 

{Ai+iBi+r) 



{A 



i+rBi) {Ai^rBi+i) ■ (Ai^rBi+r) 



(20) 



Care has to taken since the fermion creation and annihi- 
lation operators obey periodic boundary conditions when 
N/2 is odd and antiperiodic boundary conditions when 
N/2 is even. 

All that remains to do is to calculate (AiBj). We 
find the unitary transformation U that diagonalizes the 
Hamiltonian: 



4 

di 



4 = E,4^»^ 

Ci dj Uij , 



(21) 



where denotes a fermion creation operator in the old 
basis and dl denotes a creation operator in the new basis. 
Then we find that 



and in the ground state 



k<kf 



Hence 



(AB,) 
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jk- 



(22) 



(23) 



(24) 



k<kf 



This completes the expression for the longitudinal com- 
ponent of the string correlation function as a determi- 
nant. The equivalent expression for the other correlation 
functions are obtained in a very similar waycj and we 
only state the results. The transverse component of the 
string correlation function is 

i+r 



{B^A,+ l) 
(-Bi+2^i+l) 



(B^A+s) 
(Bt+2Ai+3) 



(B 



i+r-l^i+l 



A+i) {B, 



i+r-lAi+3 



(-Bi+2^i+r) 
{Bi+r-lAiJ^r) 



(25) 



Next we turn to the longitudinal component of the spin 
correlations function, which can be expressed as 



1, 



{A,B/) {AB,+r) 
{Ai+rBij (Ai^r-Bi+r) 



(26) 



the transverse component is given by 



.1, 



{B^A+l) 

{Bi+iAi+i) 



{B^A+2) 

{Bi+iAi+2) 



(Bi+r-lAi+i) (Bi^r-lAi^ 



(^r^r+.> = (4) X 

{BiA^^r) 
{Bi+lAi^r) 

■ {Bi+r-lAi+r) 



(27) 



IV. RESULTS 

In order to check the accuracy to which we can recover 
the exponents of the clean system we first display all 
correlations functions for the pure model in a log- log plot 
in Fig. 1^. In the pue« model, at zero temperature, Lieb, 
Schultz and MattiscJ showed that C^(r) is identically 
zero for r even, but decays-iwith an exponent of 2 at 
large r (odd). Later McCoycS showed that C^{r) decays 
with an exponent equal to 1/2. For the pure model 
will e qual and is the square root of (see Eq. 
(|A12| )). Hence and will decay with exponents 
1/2 and 1/4 respectively. By using a linear fit to the 
data points for the three largest system sizes the expected 
decay exponents are recovered to five significant digits 
(see Table I). 

We next add flat disorder, ^described by Eq. (^, and 
expect to see the RS exponento 2 for the spin correlations 
and C^. maps on to the transverse field Ising 
model correlatoE, which falls of with an exponent of 0.382 
in the RS phaseEJ. We therefore expect both and 
to fall off with exponent 0.382 We plot the results in Fig. 

has previously been reportedSB to show algebraic 
decay with an exponent of 2, as predicted by the RS cal- 
culation, and we confirm this result. While C'^ remains 
the same in the clean and RS case, the exponent of the 
transverse correlation function is predicted to change 
from 1/2 in the clean-case to 2 in the RS phase. Previ- 
ous numerical resultsllil for indicated a transition to 
exponential decay. Contrary to this, our results in Fig. 
^ do support an exponent of 2 in the RS case. The expo- 
nent appears to be slightly smaller than 2 for the smaller 
system sizes, but the slope gets very close to 2 for the 
larger system sizes. We believe that we do not see the 
correct exponent for smaller system sizes because has 
not yet converged to its thermodynamic limit. The last 
data point (TV — 1024) is not reliable, since for this large 
system size and strong disorder the numerical routines 
have begun to become unstable, and we do not trust the 
accuracy of this point. 
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We believe that the exponential decay seen previously 
was caused by the use of a gaussian distribution of bonds. 
As the disorder gets strong enough this will introduce 
ferromagnetic bonds into the system. The ferromagnetic 
bonds will destroy the antiferroniagnetic order leading 
to the RS phase. Westerberg et al.Q have shown that the 
spin 1 /2 AF fixed point is unstable towaifds introduction 
of FM bonds. Furthermore they showedD that the pres- 
ence of FM bonds leads to exponential decay of the spin 
correlation function. 

We note that although the decay exponent is the same 
for both components of the correlation functions, the 
prefactors differ. If the spins were coupled together in 
pairs as true singlets this would, of course, not be the 
case. We believe that the different prefactors are caused 
by residual fluctuations of the valence bonds. 

During the course of this work we did fairly extensive 
studies of various other bond distributions. In particular 
we worked with a power law distribution 



P(J) ^aJ-^+"9{J)e{l~ J). 



(28) 



The motivation for a power law distribution comes from 
the decimation RG procedure, in which strong bonds in 
the system form singlets and thereby generate weak in- 
duced couplings between their neighbors. This proce- 
dure, in effect, eliminates strong bonds and replace them 
by much weaker bonds. The energy scale is thus lowered 
and the procedure becomes exact in the low temperature 
limitldu. After renormalization the bond distribution is 
given by the above power law distribution. We believed 
that by starting as close to the fixed point distribution as 
possible it would be easier to observe the predicted RS 
behavior of the correlation functions. Another distribu- 
tion we studied was a flat distribution with mean equal 
to one, but width less than two. With neither of these 
distributions did we observe the RS behavior as clearly 
as with the distribution presented in this paper. There 
seems to be two competing factors that influence the re- 
sults. The power law distribution is very close to the 
fixed point distribution, but the bonds get distributed 
over several hundred orders of magnitude and the nu- 
merical routines become unstable and give increasingly 
unreliable results. The flat distribution with width less 
than two behaves very well numerically, but it is far from 
the fixed point, and we did not manage to reach system 
sizes large enough to observe the expected behavior. It 
appears that a flat distribution with mean one and width 
two is close enough to the fixed point for us to see the RS 
exponents, but it is not so broad that it causes problems 
for the numerical routines, except at the largest system 
size (N=1024). 

The string correlation functions decay as predicted in 
the RS phase. The decay exponent of changes from 
0.50 in the pure case to 0.382ij3 the RS phase, in agree- 
ment with earlier observationsEHI. Young and Rieger mea- 
sured the spin correlation function in the transverse-field 
Ising chain, which maps onto the transverse component 



of th e st ring correlation function in the XX model, see 
Eq. (A9). The decrease in exponent is of interest since 
the disorder actually causes the correlation function to 
fall off slower than in the clean case. We do, however, 
have to keep in mind that we are looking at the average 
correlation function and not the typical correlation func- 
tion. The average correlation function is dominated by a 
few strong couplings between distant spins. Such strong 
couplings are a characteristic of the RS phase. The typi- 
cal correlation function, on the other hand, shows a very 
different behavioiuB and will be studied below. We do 
also see the predicted change in , which changes from 
exponent 0.25 in the pure case to 0.382 in the RS phase. 

is the square of and it is interesting that the disor- 
der average of both quantities decay with the same expo- 
nent. Intuitively this can be understood by considering 
a correlation function that is dominated by a few strong 
correlations of order 1, while the rest of the correlations 
are of order 0. The actual distributions are plotted in 
Fig. I 

The typical correlation function, here obtained by ex- 
ponentiating the disorder averase of the logarithm of the 
correlation function is predictedp to decay according to 



/typ(^) « exp{~Ay/^), 



(29) 



where A is some nonuniversal constant. By plotting the 
log of C'lJyp and C'lyp against ^/r in Fig. |4|we see that this 
is indeed the case, and a linear fit to all data points gives 
the value -1.08 and -3.22 respectively for the nonuniversal 
constant A above. It is of interest to note that and 
are identical for the pure model, but in the RS phase 
and decay with different exponents. C^y-p and O^^^ 
are, however, again equivalent in the RS phase. This 
is easily understood by the mapping to the Ising chain. 
Looking at Eq. ( A12 ) it is clear that and differ in 
the RS phase since Of and Of are anti-correlated. This 
is illustrated in Fig. g and provides further evidence for 
the RS phase since in the RS phase strong bonds can 
never cross each other, and this means that if O^ is of 
order unity, then Og is bound to be small. The typical 
correlations C^yp and O^^p are equivalent since 



ln(C^(iV/2)) = ln(0;^(7V/2 - 1)) + ln(Of (iV/2 - 



1)), 
(30) 



and because if there is no dimerization in the system the 
average of In O^ {i — j) will equal the average of In Of (« — 
j). We note that since O^ = (O^)^ the logarithms of the 
typical correlations 0|yp and 0|yp will differ by a factor 
2. 

The distributions of the logarithm of the correlation 
functions 



ln(/(r))/VF 



(31) 



is supposed to sciale to a fixed distribution for large r. 
Young and Riegeill3 studied this distribution for , and 
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we present results for and in Fig. || and Fig. |^. 
As we see the distributions scale well. The distributions 
of and are very different, and it not obvious from 
looking at the plots that the typical correlations O^y-^ 
and C'^yp turn out to be the same (though they do). 

Next we focus on a dimerized system with no disorder. 
Any observable is a function of the three length scales 
in the system: the system size TV, the distance around 
the ring r, and the correlation length ^. We make the 
standard finite size scaling hypothesis and assume that 
near the critical point = oo,N = oo) we can express 
our observable in terms of a scaling function with dimcn- 
sionless arguments formed by ratios of our length scales. 
Hence a correlation function C can be expressed as 



C{N,r,O^N^f{N/r,N/0, 



(32) 



where w is a scaling exponent and f{N/r,N/^) is the 
scaling function. If we measure C for different system 
sizes, but at a fixed ratio N/r and plot CN~'^ vs. N/^ 
the curves should collapse onto a single scaling function 
f{N/^), independently of the values of N and ^. If one 
assumes that N is large enough so that C is not a function 
of N, then one can alternatively use data for one system 
size N, and different r and plot Cr~'^ versus r/f . Again 
the curves should collapse onto a single scaling function 
f{r/^), independently of the values of r and ^, but only 
as long as the value of C has reached its therraadynamic 
limit. This was done by Young and RiegeJiiJ, but in 
this paper we present data for different system sizes at 
a fixed distance half way around the ring, r = N/2, so 
that the first argument of the scaling function / is fixed 
at N/r = 2. 

Consider the string correlation function 0^{n, S) to be 
a function of the dimerization 6 and n = N/2. From 
above we know that 



0%{7i,5 = 0) cx n 4 



(33) 



and PfeutyE3 showed that (using the mapping from the 
XY model to the Ising chain studied by Pfeuty) 



0'^[n — oo,5) on Si- . 



(34) 



The correlation length for free fermions is of the form 
^ cx 5^^. Consider a string correlation function of the 
form 



(35) 



Eq. (^4|) requires that for asymptotically large values of 
X, f{x) behaves as x^, which fixes uj to —\. We thus 

expect a single curve when O^ni is plotted versus nS. In 
Fig. H a log-log plot of the scaling function shows that 
this is indeed the case. For small values of the argument 
the function is linear, f{x) oc x, but for larger values of 
the argument there is a sharp transition to the asymp- 
totic behavior, f{x) oc x*. As the dimerization gets too 
big we notice that there are corrections to scaling. 



Since the longitudinal component of the string correla- 
tion is the square of the transverse component, the above 
arguments are identical for {N, 6) , except that all ex- 
ponents equal to 1/4 are changed to 1/2. The spin cor- 
relation function also scales very well with exponent 
1/2, as shown in Fig. The longitudinal spin correla- 
tion shows fairly large corrections to the scaling (see 
Fig. |l^), but the corrections are due to the fact that 
the correlation is measured at r = A^/2 — 1 instead of at 
r = N/2. This was checked by by plotting the same re- 
sults for system sizes with N/2 odd, in which case can 
be measured at r = N/2, and there are no visible correc- 
tion to scaling. But in order to be consistent throughout 
this paper we have only shown data for system sizes with 
N/2 even. 

If we now add disorder to the dimerized system we 
expect the correlation length to be of the forma ^ c>c S~^. 
Furthermore we expect that 



0^{n,d = 0) cxn 



-0.382 



(36) 



Therefore the scaling function for both components of 
the string correlation function should read 



-0.382 



f{Nd' 



(37) 



In Fig. 



for and 



11 and Fig. |12| we plot the scaling functions 
In both cases we see a transitions from 



f{x) cx X2 to f{x) oc x"-^° . The transition is not as 
sharp as in the clean case, and it would be of interest to 
have reliable data for even larger system sizes. 

The average spin correlation functions should scale ac- 
cording to 



C=^'^(n,5) oc n-^f{nS^), 



(38) 



but since the critical correlations functions have 
not reached their asymptotic behavior at the sys- 
tem sizes considered we compensate by plotting 
C^'^(n, (5)/C^'^(n, ^ =^ vs. f{nS^), which works well, 
as can be seen in Fig. |l3| and Fig. |lj. 

The typical correlation functions are expected to scale 

as 



ftypin'^)/ftypin,S = 0) oc /(n/^typ)> 



yp^ 



(39) 



where ^typ oc ^ ^ in the RD phase. The string corre- 
lation functions and appear to follow this 
scaling behavior fairly well, see Fig. |l^ and Fig. |l^, but 
the typical spin correlation functions in Fig. y_7| and Fig. 
show quite dramatic corrections. If we assume that 
" the scaling works well, as shown in Fig. uh. 



yp 



s-1.3 



ft; 

We are not sure if this discrepancy with theory is only an 
effect of large corrections to scaling, or if there is some 
other reason. 
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V. CONCLUSION 

We have done extensive numerical simulations of the 
disordered XX model, using a mapping to free fermions. 
The decay exponents of the string and spin correlation 
functions in the RS and RD phases, predicted by RG cal- 
culation, have been found in very good agreement with 
numerical data. In particular the transverse component 
of the spin correlation function is observed to decay alge- 
braically with the correct decay expcpent, as opposed to 
exponentially as previously observedu. We explain this 
difference in terms of the bond distribution and the in- 
troduction of ferromagnetic bonds in the previous study. 
We also discuss the results of various other bond distri- 
butions, emphasizing the importance of being close to 
the fixed point distribution, but making sure that the 
numerical routines still give reliable results. Full finite 
size scaling cross over functions were presented for the 
longitudinal and transverse components of the average 
and typical spin and string correlation function in the 
RD phase, as well as for the pure dimcrizcd model. Only 
the transverse component of the string correlatioiufunc- 
tion has previously been studied in this mannerE2l. In 
appendix A various relations between the different corre- 
lation functions are derived using a mapping to the Ising 
model in a transverse field. Due to the well-known map- 
ping to free fermions, and the use of parallel computers, 
disorder averages could be taken over system sizes up to 
1024 sites, about ten times larger than in previous studies 
we are aware of. 

We are extremely grateful to Ross Hyman, Senthil To- 
dadri and Kun Yang for their many comments and gen- 
erous help. We acknowledge support from the NSF grant 
DMR 97-14055, NSF CDA-9601632 and Ella och Georg 
Ehrnrooths stiftelse. 



APPENDIX A: TRANSFORMATION TO 
DECOUPLED ISING CHAINS 

We can gain considerable insight into the behavior of, 
and relationship between various correlation functions 
by transforming the XY chain into two decoupled Ising 
chains. Following Fisher we use the following transfor- 
mation: 



rrV — QV QV 

2 2 



(Al) 



The inverse transformation is given by 

nx X X 

si= n 



(A2) 



2 2 -r-r 

qz QX ny I I y a: x 

j<n-l 

The Hamiltonian is transformed into two decoupled Ising 
chains: 



N 



N 



1=1 



(A3) 



The two chains are dual to each other and in the pure 
dimerized case one chain has coupling constant Jo and 
transverse field Je, while the other chain has coupling 
constant Je and transverse field Jq. 

Next we look at how the string correlation functions 
transforms: 

0-(r) = {S^S:^, ■ ■ ■ Sr+.) = {aW^^^r^- (A4) 

If r is even then crf_x and crf^^ will belong to differ- 
ent chains and the expectation value < af_iaf^j, >=< 
af_i >< af_^_^ > will vanish. If r is odd 0^{r) trans- 
forms to the Ising correlator: 0^{r) — {a'f_iaf_^_^) . The 
y component becomes (for r odd): 

Oy{r) = {SfSf^, ■ ■ ■ = (af • • • ^f+r-i)- (A5) 

The expectation values of and have to be identical 
and hence 



(A6) 



This relates the expectation value of a correlation func- 
tion on one chain to the expectation value of another 
correlation function on the other chain. 

The z component of the string correlation function 
transforms (r odd) as 



' i+r-l" i+rl 



(A7) 



Using the above equality we find that 

= {<jU<y^+r? = 0^r)\ (A8) 

and interestingly enough the transverse component of the 
string order is simply the square of the longitudinal com- 
ponent. 

Turning to the spin correlation function we find that 



C%r) = = (a,_ia,a,+,_ia,+,). (A9) 



Hence 



C^(r) = 



< (T,,_icri+r >< o'iCTj+j.-i >, if r odd 

< (Tj-iCTi+r-i >< CTiCTi+r >, if r even. 



j<n-l 



(AlO) 
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The criterion that C = gives us back Eq. (|Ag). The 
z component of the spin correlation transforms to 



1 "i±2"i+4 



' i+r 



> X 



> X 



< erf fff ■ 



if r odd 



if r even 



(All) 



To summarize we want to emphasize a few important 
relations between various correlation functions. These 
equalities follow directly from the fact that the XY chain 
separates into two decoupled Ising chains. The expres- 
sions are obtained by using the above transformations 
and assuming that N/2 is even, which is the case in our 
simulations. The relations are: 

C^(iV/2) = Ol{N/2 - l)Of (iV/2 - 1) 
0:{N/2-l)^iO:iN/2^l))^ (A12) 
0:(A^/2-l) = (0:(A^/2-l))2. 
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FIG. 1. The spin and string correlations functions in the 
clean XX model. System sizes go from = 32 to = 2048 
and n = N/2. The solid lines are linear fits to the three 
largest systems. Slopes are given in Table |. 
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FIG. 3. The distribution of the transverse and longitudinal 
components of the string correlation function in the RS phase 
for system size A'^ = 256. 
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FIG. 4. The typical correlation functions plotted versus the 
FIG. 2. The spin and string correlation functions in the square root of the system size. The solid lines are linear fits 

disordered XX model. Disorder averages are taken over to all data points. 

5 X 10^ — 10^ configurations, and system sizes vary from 
= 32 to = 1024. The solid and dashed lines have slopes 

-2 and -0.382 respectively. 
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FIG. 13. The longitudinal component of the spin correlar 
tion function in the dimerized and disordered XX model. 
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FIG. 16. The transverse component of the typical string 
correlation function in the dimerized and disordered XX 
model. 
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FIG. 14. The transverse component of the spin correlation 
function in the dimerized and disordered XX model. 




FIG. 17. The longitudinal component of the typical spin 
correlation function in the dimerized and disordered XX 
model. 
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FIG. 18. The transverse component of the typical spin cor- 
relation function in the dimerized and disordered XX model. 




FIG. 19. The transverse component of the typical spin cor- 
relation function in the dimerized and disordered XX model, 
with ^ oc 5^-^. 
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